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ABSTRACT 

Context. Asteroseismic determinations of structural parameters of hot B subdwarfs (sdB) have been carried out for more than a 
decade now. These analyses rely on stellar models whose reliability for the required task needs to be evaluated critically. 
Aims. We present new models of the so-called third generation (3G) dedicated to the asteroseismology of sdB stars, in particular to 
long-period pulsators observed from space. These parameterized models are complete static structures suitable for analyzing both p- 
and ^-mode pulsators, contrary to the former second generation (2G) models that were limited to p-modes. While the reliability of 
the 2G models has been successfully verified in the past, this important test still has to be conducted on the 3G structures. 
Methods. The close eclipsing binary PG 1336—018 provides a unique opportunity to test the reliability of hot B subdwarf models. We 
compared the structural parameters of the sdB component in PG 1336—018 obtained from asteroseismology based on the 3G models, 
with those derived independently from the modeling of the reflection/irradiation effect and the eclipses observed in the light curve. 
Results. The stellar parameters inferred from asteroseismology using the 3G models are found to be remarkably consistent with both 
the preferred orbital solution obtained from the binary light curve modeling and the updated spectroscopic estimates for the surface 
gravity of the star The seismology gives = 0.47 1 ± 0.006 Mg , = 0. 1474 ± 0.0009 R,^; , and log g = 5 .775 ± 0.007, while orbital 
modeling leads toM* = 0.466±0.006 M.y, = 0.15 ±0.01 Rq,, logg = 5.77 ±0.06, and spectroscopy yields logg = 5.771 ±0.015. 
In comparison, seismology from a former analysis based on the 2G models gave very similar results with = 0.459 ± 0.005 Mq, 
= 0.151 ±0.001 Rr.^, and logg = 5.739 ±0.002. We also show that the uncertainties on the input physics included in stellar models 
have no noticeable impact, at the current level of accuracy, on the structural parameters derived by asteroseismology. 
Conclusions. The stellar models (both of second and third generation) presently used to carry out quantitative seismic analyses of 
sdB stars are reliable for the task. The stellar parameters inferred by this technique, at least for those that could be tested (M*, R, and 
logg), appear to be both very precise and accurate, as no significant systematic effect has been found. 

Key words, stars: subdwarfs - stars: oscillations - stars:interiors - stars: binaries: eclipsing - stars: individual: PG 1336—018 

1. Introduction Major uncertainties remain about the details of the forma- 
tion of the hot B subdwarfs. One important clue is the fact that 

Hot B subdwarf (sdB) stars are hot and compact objects with the fraction of sdB stars in binai'ies with stellar c ompanion is 

^ effective temperatures T^s between 20 000 - 40 000 K and sur- found to be relatively high, around - 50% (see, e.g., AilaidetaD 

Pi . face gravities log g in the range 5.0 - 6.2 (see, e.g.. Green et al. 1 994; Maxted etal...2001; Green et al 2008). Plausible forma- 

I2OO8). They occupy the so-called Exti-eme Horizontal Branch tion channel s were modeled in det ail from binary population 

(EHB), burning helium in the core and having a very thin resid- synthesis by i Han et all (l2002l l2003h and include binary evolu- 

ual hydrogen-rich envelope. The sdB stars dominate the pop- tion via a common envelope (CE), stable Roche lobe overflow 

ulation of faint blue objects down to V - 16 and are mainly (RLOF), and the merger of two helium white dwarfs, resulting 

found in the thick galactic disk, but also in the halo, as glob- in a single sdB stai-. The proposed scenario for forming single 

ular clusters members a nd in the galactic bulge (|Altmann et alj stars remains unsettled however, as fliey could also result from 

|2004|Busso et al.||2005). Their presence is indirecdy revealed in single star evolution through enhanced mass los s at the tip of 

distant objects flirough their copious UV emission that is respon- the RGB at the moment of He-burning ignition dPorman et al.1 

sible for the UV-uptum phenomenon (an excess of UV emission) (1993 ). or through a delayed He-flash, at higher effective tem- 

obser ved, for instance, in old elliptical galaxies (Brown et all peratures during the collapsing phase after having left the RGB 

|1997|). (iD'Cruz et al.ll996h . Another possibility is currently gaining ob- 
servational support and involves the interaction of the star with 

Send offprint requests to-.W.W^n Grootel substeUar companions (brown dwarfs and/or planets) in 
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a common envelope evolut ion du r ing th e red giant phase. This 
scenario first proposed by 'SokerJ (Il998l) is strongly suggested 
by the discovery of small, nearly Earth-sized planets around the 
isolated sdB star KPD 1943+4058 (KOI 55) that could be rem- 
nants of former giant planets that survived engulfment in the 
red giant envelope and contributed to expel the stellar envelope 
(ICharpinet et alJ 1201 lai) . These distinct evolutionary scenarios 
leave a clear imprint on the binary fraction and distribution of 
sdB and their companion stars (CE evolution produces sdBs in 
very close binary systems preferably with white dwarfs com- 
panions, RLOF gives rise to longer period sdB + main-sequence 
star binaries), but also, most importantly, on the stellar mass 
distribution for sdB stars. After core-He exhaustion, the sdB 
stars evolve directly t oward the white dw arf cooUng sequence 
as AGB-Manque stars dPorman et al 

Asteroseismology should help clarify the question of the for- 
mation of sdB stars. The sdB stars indeed host two groups of 
nonradial pulsators. The first group, named the V36 1 Hya (orig- 
inally EC 14026) stars dKilkennv et al.lll997h . exhibits short- 
period pulsations in the range 80 — 600 s with amplitudes of 
a few milli-magnitudes. The pulsations correspond to low-order 
low-degree p-modes, which, in sdB sta rs, have significant a mpli- 
tudes mainly in the outermost layers ( Charpinet et al.l2000l) . The 
second group, discovered by iGreen et al. (20()3h and referred 
to as VI 093 Her stars, exhibits much longer periods from 45 
min to a few hours. This corresponds to mid- and high-order 
low-degree g-modes, which can propagate in deep regions of 
the star, down to the convective He-burning core. Both p- and 
^-mode pulsations are driven by the same fc-mechanism pow- 
ered by local envelope accumulat ions of heavy elements such 
as iron due to radiative levitation (Charpinet et al."l996l Il997t 
[Fontaine et al. 2003; see also Hu et al. 2011). Quantitative aster- 
oseismology of pulsating sdB stars has developed considerably 
during the past ten years, giving us new ways of revealing the 
properties of these objects at unprecedented level of accuracy. 
This includes, among many other structural parameters, tight 
constraints on their masses. The re s ults ob tained so far are com- 
piled in Table 1 of iFontaine et al.l (1201 2|) . with the appropriate 
reference to each analysis having been carried out. 

The method has long relied on stellar envelope models 
of the so-called second generation (2G) that incorporate the 
nonuniform abundance profiles of iron predicted by equilib- 
rium between radiativ e levitation and gravitational settling (see 
ICharpinet et al.ll2002l for details). These models are suitable for 
asteroseismology of pure p-mode pulsators whose oscillation 
modes are sensitive only to the details of the outermost regions 
of the star A very important test of accuracy for the seismic 
results obtained on the basis of these models w as made using 
the star PG 1336-018 (ICharpinet et al.l l2008bl) . This object is 
one of the only two known sdB + dwarf M close eclipsing bi- 
naries where the sdB component is a pulsating sdB star (the 
other one being 2M 1938+4603 recently discovered with Kepler, 
I0stensen et al , .2010) . PG 1336-018 permitted to compare re- 
sults obtained from the two independent techniques of astero- 
seismology and orbital light curve and eclipses modeling. It was 
found that asteroseismology is remarkably c onsiste nt with one of 
the orbital solutions uncovered by Vucko vic et al.l (i2007h . This 
close agreement showed that the seismic inference is both pre- 
cise and accurate (i.e., free of significant systematics), at least for 
the stellar parameters that could be tested. This was of course an 
important finding attesting that the 2G models used for the quan- 
titative asteroseismic studies of V361 Hya stars are sufficiently 
reliable for the task. 



Recently, the advent of the space missions CoRoT and 
Kepler opened up new opportunities to gather sufficiently high 
quality data for the very low-amplitude, long-period g-mode pul- 
sations seen in sdB stars for asteroseismological purposes. Three 
^-mode sdB pulsators have indeed already been analyzed on this 
basis, bringing for the first time information on their He-burning 
cores - including size and composition - along with the determi- 
natio n of other structur al parameters (Van Grootel et al. 20 lOa^^ 
Cha rpinet et al.ll201 Ibl) . Because the 2G envelope models are no 
longer suitable for quantitative analyses of ^-mode pulsators, a 
third generation (3G) of stellar models for accurate evaluations 
of the g-mode pulsation periods had to be developed for that pur- 
pose and were used in the studies cited above. These are com- 
plete stellar structures (as opposed to the 2G envelope models) 
that include a detailed description of the central regions, as well 
as, like their 2G counterparts, the nonuniform abundance profile 
of iron in the envelope predicted by microscopic diffusion. 

While the 2G models passed the test of reliability for astero- 
seismology of p-mode sdB pulsators with the case of PG 1336- 
018, the new 3G models still have to be verified in a similar 
way. In the following section (Section 2), we discuss the need for 
parameterized models in quantitative asteroseismology, and we 
describe the new 3G models that we use for analyzing pulsation 
data obtained on hot B subdwarfs. In Section 3, we reanalyze the 
available data on PG 1336-018 with these 3G models in order 
to confront the seismic inferences obtained on the fundamental 
parameters of this star with those resulting from 1) the former 
seismic analysis based on 2 G models, and 2) the independent 
orbital light curve analysis of lVuckovic et al.l (l2007l) . This test is 
important for linking older asteroseismic results obtained with 
2G models with newer results based on 3G models in order to 
estimate potential systematics. We also carry out experiments in 
order to assess the importance of known model uncertainties on 
the stellar parameters determined by asteroseismology. A sum- 
mary and conclusion are finally provided in Section 4. 

2. Third generation of sdB models for seismology 

2.1. Background 

It may be worthwhile here to give a brief historical account of 
how our group came to the development of successive genera- 
tions of models for the specific purpose of applying asteroseis- 
mological techniques to sdB stars. We refer to as first gener- 
ation (IG) models those stellar equilibrium structures for sdB 
stars with a uniform metallicity. They include full evo lutionary 
models such as those provided by Dor man et al.l (1 19931) at a time 
when no pulsating sdB stars were known, and when one of us 
(S.C.) had started a Ph.D. thes is to investigate the asteroseismo- 
logical potential of these stars dCharpinetl 1 9991) . They include as 
well static structures made of an envelope (calibrated in part on 
the basis of the evolutionary models of Dorman) sitting on top 
of an inert ball. Those envelope models allowed quite a bit of 
flexibility, in particular for investigating the effects of changing 
the metallicity on the overall stability of pulsation modes. It was 
thus found that uniform metallicity models (both evolutionary 
and static ones) are generally stable against pulsations, unless 
the metalhcity is incre ased to some unreaUstically high values 
(ICharpinet et al.lll996h . 

These investigations with IG models identified a region of 
local instability associated with an opacity feature (the so-called 
Z-bump) even though the modes were globally stable. Given that 
sdB stars are very hot objects, and given that they are known 
to be chemically peculiar, radiative levitation was the obvious 
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new ingredient to incorporate into a next generation of equilib- 
rium structures. The idea was to verify if enough iron could ac- 
cumulate in the driving region through radiative levitation and 
push some pulsation modes toward instability. Hence, we de- 
veloped our second generation (2G) models, which are again 
static envelope structures, but which incorporate a nonuniform 
distribution of iron (a representative Z-bump element) as deter- 
mined by microscopic diffusion. In detail, the iron profiles are 
obtained through detailed calculations of radiative forces under 
the assumption of an equilibrium between radiative levitation 
and gravitational settling (resulting from the combined action 
of the pressure gradient and electric fields). This is called here- 
after "diffusive equilibrium", disregarding, at this stage, other 
competing processes (see Section 13.5.1b . The computations of 
the radiative accelerations themselves and of the equilibrium 
ab undances of Fe are done following the procedure described 
in Chav er et al.l (11995*). These 2G models have been extremely 
successful at explaining the location of the p-mode sdB pul- 
sators in the log g-T^ff plane as well as their nonadiabatic prop- 
erties, and they have remarkably well resisted the test of time 
dCharpinet et al. 1997, 2001, 2006; Randall et al. 2012). The 
building assumption that the local abundance of Fe in a sdB 
star could be estimated by equating the effects of radiative lev- 
itation with those of gravitational settling was ultimately con- 
firmed through detailed time-dependent calculations (see, e.g., 
[Fontaine et al. 2006; Charpinet et al. 2008ai). 

Our driving motivation in the early years of sdB seismology 
was to explain the very existence of pulsators and account for the 
observed ranges of excited periods in these stars. Hence, "sim- 
ple" envelope structures (such as the 2G models) incorporating 
the appropriate details of the driving/damping zone proved to 
be perfectly suitable for investigating the non adiabatic proper- 
ties of pulsating sdB stars. At the same time, iCharpinetl (1 19991) 
ca rried out a very detailed analysis with the evolutionary models 
of iDorman et al.l (fl993), and came to the conclusion that sim- 
ple envelope models (obtained by cutting off the central zone 
of Dorman's evolutionary models or by building them with the 
same constitutive physics) would also be quite useful for quan- 
titative adiabatic asteroseismology as long as only p-modes are 
involved dCharpinet et al.ll2002ll) . In contrast, systematic effects 
(the periods are increasingly overestimated with increasing ra- 
dial order, reaching up to differences of 10—20%) appear when 
dealing with g-modes, which is not surprising given that the cen- 
tral region of a sdB model carries a significant weight in terms 
of period formation for those modes contrary to p-modes. The 
use of the 2G models thus allowed us to combine both adi- 
abatic and nonadiabatic seismology for quantitatively investi- 
gating the properties of the /9-mode pulsators. Our initial ef- 
fort in that direction was the study of the star PG 0014H-067 
which led to the first successful seismic inferences on the global 
structural parameters of a sdB pulsator, including its total mass 
dBrassard et al.l 1200 Ih . This was before the discovery of the 
cooler and lower gravity long-period g-mode pulsators in 2003. 
The need for complete stellar models (incorporating the impor- 
tant central region for g-modes) became apparent when suffi- 
ciently high quality data on g-mode pulsators became available 
from space missions. 

2.2. Description 

The 2G models are envelope structures built under the assump- 
tion that the luminosity is constant throughout the entire enve- 
lope. This, it turns out, is an excellent approximation for He- 
core burning stars on the EHB, as can be seen in Fig. 2 of 



ICharpinet et al.l (l2000h . for example. As a consequence, only 
three equations remain to define the structural properties of such 
an envelope, and the integration is carried out from the sur- 
face inward down, typically, to a fractional mass depth log^ = 
log(l — M(r)/M*) — —0.1. There are four primary parameters 
to specify a 2G model: the effective temperature T^g, the sur- 
face gravity log g, the total mass of the star M*, and the mass 
contained in the H-rich envelope, usually expressed as a loga- 
rithmic fractional mass, log(Menv/A^*)- The first two (T^ff and 
log g) can be directly related to observable quantities obtained 
by spectroscopy. 

In the case of the 3G models, these are now complete stel- 
lar structures assumed to be in strict thermal and mechanical 
equilibrium. Thus, the luminosity of such a static model is ex- 
actly balanced by core He burning (and the very small contri- 
bution of H burning at the base of the H-rich envelope). The 
four usual equations of ID stellar structure (minus the time- 
dependent term) are simultaneously solved with the help of a 
relaxation scheme. The physical conditions at the outer bound- 
ary are handled through a gray atmosphere scheme. A thin con- 
vection zone is usually created around the opacity peak in the en- 
velope of a sdB model, but the flux carried by convection across 
that thin zone is quite small, and well below the base of the at- 
mosphere (Fig.[Tli. When present, convection is handled through 
one or another of the various versions of the mixing-length the- 
ory. In the present application, convection is modele d in terms 
of the standard MLl version of lBohm-Vitensd dl958l) . 

The constitutive physics incorporates the nuclear reaction 
rates of|Caughlan & Fowler (1988) for H and, particularly. He 
b urning. Electron screening is handled through the prescription 
of lGraboske et al.l(ll973i) . Since it was developed to model white 
dwarf stars as w ell, it also incorpo rates neutrino production rates 
as described by lltoh et all (Il996l) . and references therein, and 
conductive opacities, also made available by Itoh et al. (200^, 
and references therein. The radiative opacities are OPAL data of 
1995, including additional tables for iron in a hydrogen back- 
ground especially computed for our needs by Forre st Rogers in 
1996 and briefly descrilDed in lCharpinet et al.l (1 19971) . These spe- 
cial tables are essential for properly mapping the opacity profile 
of Fe associated with radiative levitation. As indicated above, the 
nonuniform distribution of Fe with depth is obtained under the 
assumption of a perfect equilibrium between radiative levitation 
and gravitational settling following the procedure described in 
Chaver et al. ( 1995). The equation of state (EOS) is an in-house 
construction covering a large domain of the temperature-density 
plane, sufficiently large to cover all types of stars, including 
white dwarfs. Three regimes are involved. For the low-density 
regions, a network of Saha equations is solved for a mixture of 
radiation and an almost ideal (including Coulomb corrections), 
nondegenerate, partially ionized gas composed of a mixture of 
H, He, C, and O in various proportions as needed. In the region 
of partial ionization where noni deal and dege neracy effects are 
important, we use the EOS of iSaumon et al. (1995 ) for H and 
He, an improved version of the EOS of Fontaine et al.l (fT977l) 
for C, and an unpublished EOS for O developed at Universite de 
Montreal in 1998. The third regime covers the high-density do- 
main corresponding to the fully ionized plasma in the liquid and, 
ultimately, in the sohd p hase according to the basic physics de - 
scribed bv lLambI (119741) and improved by Kitsikis et all d2005l) . 
Efforts were made to smoothly bridge the different EOS sur- 
faces across the boundaries of the different regimes. This is of 
particular importance for pulsation studies. Let us note, how- 
ever, that the third regime covering the high-density domain 
is not relevant here for sdB stars. Interpolation in composi- 
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Fig. 1. Some structural properties for a representative 3G model. 
Upper panel: run of the logarithm of the square of the Brunt- 
Vaisala frequency as a function of fractional mass depth (solid 
curve) and equilibrium abundance profile of Fe (dotted curve, 
shifted upward by 1.5 dex). The solar ratio is log Fe/H = —4.5. 
(—3.0 on the shifted scale). Lower panel: run of the logarithm 
of the opacity as a function of depth (solid curve), and profile of 
the He distribution (dotted curve). In both panels the locations 
of the atmospheric layers, defined as those with an optical depth 
between log T = 1.0 and log T = —1.0 are indicated by dots. 

tion is handled follo wing the additive volume prescription of 
iFontaine et"al] ( ll977h . 

The primary parameters to define a 3G model are: the total 
mass of the star M*, the mass contained in the outer H-rich enve- 
lope log(Menv/A^*), the mass contained in the (usually convec- 
tive) core log(Mcore/A^*), and the chemical composition in the 
core (with the constraint X(He) + X(C) + X(0) + Z = 1.0). The 
stratified overall structure consists of a core, surrounded by a He- 
rich mantle, itself surrounded by a H-rich envelope in which Fe 
radiatively levitates under the assumption of diffusive equilib- 
rium. The composition transition layers at the interface between 
the H-rich envelope and the He mantle, and those at the interface 
between the He-rich mantle and the core are modeled in terms 
of two additional parameters which become of relevance in pres- 
ence of diffusion. In the present application, the composition 
transition zones have been calibrated on the basis of full evolu- 



tionary models such as those of iDorman et al ] (119931) . which do 
not take diffusion or possible sources of turbulence into account. 
These extra parameters are therefore held fixed in the present ex- 
periment (see, however, Section l3.5.2b . Likewise, the value of Z 
is fixed at 0. 

As an illustrative example. Fig. [T]depicts some of the struc- 
tural properties of the 3G model identified below as the model 
that best accounts for the available seismic data on the pulsat- 
ing sdB component of the PG 1336—018 system. The relatively 
sharp peaks which characterize the profile of the Brunt- Vaisala 
frequency near log q ~ —0.35 and log q ~ —3.8 are clearly as- 
sociated with the abrupt change of chemical composition in the 
composition transition zones. The signature of these transition 
zones is also evident on the opacity profile, although the mag- 
nitude of that signature is much less important in that case. The 
model features a convective core and a much less significant con- 
vection zone in the outer envelope, in the sense that the latter 
carries at the most only '^0.6% of the flux. This very weak con- 
vection zone is clearly associated with the maximum of the opac- 
ity profile, but does not change in any significant way the basic 
character of the driving mechanism, which remains essentially a 
pure fc-mechanism. It is the location and shape of the Fe opacity 
peak that provide the conditions for exciting p-mode pulsations 
in the model. Again, radiative levitation is the key element for 
boosting the abundance of Fe in the driving region, which is nec- 
essary for destabilizing pulsation modes. The opacity profile in 
the envelope not only affects the driving/damping processes, but 
also the adiabatic properties (particularly the pulsation periods 
for /?-modes) through its direct effects on the mechanical struc- 
ture of the model. Finally, this 3G model is a complete model 
and, therefore, provides a reliable equilibrium structure for de- 
riving not only p-mode periods, but also g-mode periods. 

2.3. The need for parameterized models 

When we started investigating the application of the forward as- 
teroseismological method to pulsating hot subdwarf and wh ite 
dwarf stars a decade ago (see, e.g., iBrassard et al.l 12001 1 or 
Fontaine & Brassard 2002), we realized at the outset that the use 
of full evolutionary models would not be practical. This is be- 
cause it is nearly impossible, even with large computer clusters, 
to cover finely all of the relevant domains of parameter space 
with evolutionary sequences; it would take too much computing 
time. The practical consequence is that evolutionary sequences, 
necessarily limited in number, may actually miss the correct re- 
gion of parameter space where resides the best seismic model. 

As a fallback position, we turned to less realistic, static mod- 
els of the kind described here. Those allow maximum flexibility 
for very exhaustive searches in parameter space and the eventual 
identification of an optimal seismic model as the case may be. 
Of course, such a seismic model must ultimately be validated by 
detailed evolutionary calculations. The parameterized model ap- 
proach can also be useful for identifying readily the dependences 
of the pulsation periods on specific model parameters. They can 
also help identifying rather directly the parts of the input consti- 
tutive physics that may need improvement. 

The parameterized static 3G models that we discuss in this 
paper are still crude models of hot subdwarf B stars. We be- 
lieve, however, that they are defined in terms of the most sensi- 
tive parameters from a seismic point of vie w. It is important that 
they be tested as thoroughly as possible. Brassard & F ontain3 
(2008) presented the results of a first test connecting the seismic 
mod el of PG 0014 - H067 derived from 2G structures as described 
by .Brassard et al.l (120011) with the new seismic solution based 
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on the first 3G models that were being developed for that star. 
Interestingly and encouragingly, both solutions turned out to be 
essentially the same. A still more stringent test is described in 
the next section. 



3. Seismic test of accuracy for the new models 

3.1. PG 1336-018: A Rosetta stone for asteroseismology 

The system PG 1336—018 (NY Virginis) is one of the very few 
known sdB + dM close eclipsing binaries, and one of the only 
two known objects of this type where the sdB component is a 
pulsating star ( Kilkenny et al. 1998). This makes PG 1336-018 
the equivalent of a Rosetta stone permitting independent tech- 
niques used to derive some of the stellar parameters (namely, 
orbital modeling, asteroseismology, and spectroscopy) to be di- 
rectly confronted. One can therefore check the reliability of the 
methods and models, as well as the presence (or lack thereof) 
of biases in the derived parameters. This rare opportunity pro- 
vides extremely important information on the accuracy actually 
achieved, as generally only the internal precision of the measure- 
ments given by the error estimates can be obtained. 

Parameters for the two components of PG 1336-018 have 
been determined by combining multicolor ULTRACAM/VLT 
light curves and the radial velocity (RV) curve obtained from 
UVESA'^LT spectra (Vuckovic et al. 2007). These authors mod- 
eled the reflection/irradiation effect and eclipses seen in the 
light curves and found three solutions - each deriving among 
other parameters a mass, a radius, and a surface gravity for the 
sdB star - of equal statistical significance, due to the large pa- 
rameter space and correlations between some parameters (see 
Table |2] for a summary of the relevant values). Independently, 
the fir st asteroseismic analysi s of the sdB component carried 
out bv lCharpinet et al. (l2008bh was base d on the pulsation spec- 
trum provided bv Kil kennv et al] (l2003h . These authors identi- 
fied 28 periodicities in the 96 — 205 s range by analyzing the 
multisite Whole Earth Telescope (WET) campaign held in April 
1999, which constitutes the most extensive photometric monitor- 
ing available to date for this star The seismic analysis was based 
on the 2G envelope models and yielded a soluti on with param- 
eters f alling very close to one of the solution of IVuckovic et alJ 
(l2007h . The latter derives for the sdB component a total stellar 
mass of = 0.466 ± 0.006 M© and a radius of R, = 0.15 ± 
0.01 Rq, while the seismic analysis of the observed p-modes led 
to a total stellar mass of M* = 0.459 ± 0.005 Mq and a radius of 
/?* = 0.151 ± 0.001 Rq. This close agreement is within the lc7 
error estimates of each technique. 

In the following, we present a new asteroseismic analysis 
of the sdB component of PG 1336-018 based on the same set 
of observed frequencies, but using, this time, the 3G complete 
models previously described. We also take this opportunity to 



update the spectroscopic estimat es of the atmospheric p arame- 
ters of PG 1336-018 reported in Charpinet et al.' (12008b') on the 
basis of better spectra and of new NLTE model atmospheres in- 
cluding metals (subsection [321 see also TableflJ. Our goal is to 
provide a thorough comparison of the various results and test the 
validity of the 3G models for asteroseismology of sdB stars. 



3.2. Improved spectroscopy for PG 1336-018 

There are two reasons why we were prompted to reinvestigate 
this issue. First, in view of the high importance of PG 1336—018 
as a testbed for stellar modeling and asteroseismology, we kept 
obtaining spectra over the orbital phase and ended up with a 



combined zero-velocity optical spectrum reaching a very high 
val ue of S/N ~ 491 around 4000 A (compared to 175 in 
iCharpinet et al.ll200 8b). We refer to that spectrum as Bok9 in 
what follows, reflecting the fact that it was gathered at the 2.3- 
m Bok Telescope of the Steward Observatory Kitt Peak Station, 
and that it has a resolution of 8.7 A at 4000 aQ. It is the sum 
of 14 individual spectra, 6 obtained when the reflection effect 
in the system was near minimum, and 8 taken at other orbital 
phases. We also combined the six former spectra into a single 
one named Bok9L (as in "low" reflection effect; S/N ^ 344 in 
the blue), and the eight others to obtain Bok9H (S/N 350 in the 
blue) in order to investigate the effects, if any, of the irradiation 
of the secondary on the derived atmospheric parameters. 

The second reason is that we have now developed the capac- 
ity to build NLTE model atmospheres and synthetic spectra with 
arbitrary heavy element abundances for a wide range of H and 
He compositions within reasonable computation times. It is well 
known that hot B subdwarfs are all chemically peculiar stars and 
that their atmospheric compositions are nonstandard and vary 
from star to star The effects of metals on the determination of 
the atmos pheric pararneters of sdB stars have been studied in the 
past (e.g.. lHeber et aljfeOOOl) . but this has been done only within 
the LTE approximation. In the present case, we use two of our 
recent NLTE model grids to investigate potential effects in the 
context of PG 1336-018. 

The firs t grid is similar to the NL TE, H and He model spec- 
tra used in ICharpinet et al.l (l2008bl) . except that our new syn- 
thetic spectra extend out to the red to include the Ha region 
(not previously considered). The second NLTE grid includes a 
representative heavy el ement mixture inspired from the work of 
iBlanchette et al.l (l2008h . These authors presented a NLTE anal- 
ysis of FUSE spectra to derive the abundances of several astro- 
physically important elements in the atmospheres of five typi- 
cal long-period pulsating sdB stars. The five stars analyzed by 
IBlanchette et al.l (l2008l) show very similar abundance patterns, 
and it is from these results that we derived a representative com- 
position using the most abundant heavy elements. Hence, we as- 
sumed atmospheres containing C (1/10 solar), N (solar), O (1/10 
solar). Si (1/10 solar), S (solar), and Fe (solar). These two 3D 
grids (reff, log g, and log A^(He)/A^(H)) were computed with the 
pubhc codes TL USTY and SYNSPEC (iHubenv & Lanzlll995l: 
iLanz & Hubenvlll995l:lLanz et alJll997^ . These codes have been 
parallelized by one of us (PB.) to efficiently run on a computer 
cluster With these two grids of synthetic spectra, as well as a 
corresponding LTE grid for comparison purposes, we analyzed 
the spectra Bok9, Bok9L, and Bok9H. We also reanalyzed the 
blue spectrum PB6, renamed Bok6 from then on, as it has a 
higher resolution (^ 6 A), but a somewhat lower S/N of ^ 80. 

We summarize, in Table [T] the results of several analyses 
of the Bok6 and Bok9 spectra, including a test with a metal- 
free LTE grid for comparison purposes. For practical use (espe- 
cially in the seismic investigation that follows), the solutions ob- 
tained with the NLTE grid including metals are to be preferred. 
Figures |2] and [3] show the spectral fit obtained with that grid for 
the Bok6 and Bok9 spectra, respectively. Given the exceptional 
S/N ratio of that latter spectrum. Fig. |3]reveals an extraordinary 
good match, notwithstanding the He II 4686 A line (although we 
had to increase the vertical scale by a factor of 5 in that case or, 
otherwise, the differences would not have been seen in the plot). 
Except for the redder regions where the S/N has significantly 
decreased (specifically, near the Ha and He I 6678 A lines), the 



' The spectrum named BG9 in lCharpinet et al.l ( l2008bl) represents the 
sum of our first few spectra and is, therefore, part of Bok9. 
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Table 1. Summary of all available atmospheric parameters of PG 1336—018 derived from spectroscopy. 



Spectrum 


Type of model 


Teff (K) 




log A'(He)/A'(H) 


Reference 


SAAOl 


LTE, H 


33,139 ± 1000 


5.78 ±0.10 




Kilkenny et al. (1998) 


SAA02 


LTE, H 


32,895 ± 1000 


5.67 ±0.10 




Kilkenny etaL (1998) 


VLT/UVES 


LTE, H, He 


31,300 ±250 


5.60 ± 0.05 


-2.93 ± 0.05 


Vuckovic et al. (2007) 


PB6 


NLTE, H, He 


33,220 ± 170 


5.75 ± 0.04 


-3.11 ± 0.25 


Charpinet et al. (2008) 


BG9 


NLTE, H, He 


32,380 ± 150 


5.77 ± 0.03 


-2.87 ± 0.17 


Charpinet et al. (2008) 


Bok6 


NLTE, H, He 


32,874 ± 148 


5.705 ± 0.028 


-3.066 ± 0.250 


This work 


Bok6 


NLTE, H, He, metals 


33,285 ± 169 


5.707 ± 0.028 


-3.016 ± 0.216 


This work 


Bok9 


LTE, H, He 


31,959 ±71 


5.823 ± 0.018 


-2.847 ± 0.101 


This work 


Bok9 


NLTE, H, He 


32,359 ± 84 


5.793 ±0.017 


-2.879 ± 0.095 


This work 


Bok9 


NLTE, H, He, metals 


32,649 ± 97 


5.794 ± 0.017 


-2.892 ± 0.098 


This work 


Bok9H 


NLTE, H, He, metals 


32,608 ± 104 


5.798 ±0.018 


-2.839 ± 0.098 


This work 


Bok9L 


NLTE, H, He, metals 


32,748 ±111 


5.800 ± 0.020 


-2.893 ± 0.113 


This work 


Former adopted value 


NLTE,H,He 


32,780 ± 200 


5.76 ± 0.03 


-2.94 ±0.14 


Charpinet et al. (2008) 


New adopted value 


NLTE,H,He,metals 


32,807 ± 82 


5.771 ±0.015 


-2.918 ± 0.089 


This work 




Fig. 2. Model fits (thick curves) to the hydrogen and helium lines (thin curves) available in our time-averaged, mid-resolution 
spectrum Bok6. 



other weak features seen in the observed spectrum are real. For 
instance, the weak structures in the red wings of H7, H5, and 
He are due to metallic absorption of unknown origin. Likewise, 
the absorption structure seen in the red wing of He I 5876 A is 
also real and is due, in that case, to weak interstellar reddening 
associated with the well known Na I doublet. For our purposes, 
we will adopt the weighted means of the atmospheric parame- 
ters as inferred in Fig. |2]and|3] Hence our updated estimates of 
these parameters for the sdB component of PG 1336—018 are 
Teff = 32 807 ± 82 K, log g = 5.771 ± 0.015 (in cgs units, as 
elsewhere in the text and in the tables), and log A^(He)/A^(H) = 



—2.918 ± 0.089. These estimate s are ve ry close to the former 
values adopted in lCharpinet et al.l (l2008bl) . but are more precise. 

Finally, we note, from Table [T] that the solutions obtained 
with Bok9H and Bok9L are, well within the uncertainties, the 
same as those obtained with the full spectrum using the same 
grid of models. This shows that the inferred atmospheric param- 
eters do not depend in any significant way on the orbital phase 
at which the spectra are taken, i.e., the results are not affected by 
the irradiation and reflection effects on the secondary. Moreover, 
the NLTE with metals vs NLTE without metals comparison for 
both the Bok6 and Bok9 spectra indicates that the presence of 
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Fig. 3. Model fits (thick curves) to the hydrogen and helium lines (thin curves) available in our time-averaged, very high signal-to- 
noise, low-resolution spectrum Bok9. 



metals (at least at the level that we considered) does not strongly 
change the values of the inferred atmospheric parameters. We 
find only a small systematic effect in effective temperature (the 
inclusion of metals slightly increases it), and negligible effects 
for the two other quantities. This shows that our spectroscopic 
determination of the atmospheric parameters of PG 1336-018 is 
very robust. 

3.3. Parameters from asteroseismology using 3G models 

The procedure to derive the structural parameters of the pulsat- 
rn g sdB star in PG1336-0 18 is the same as the one described 
in ICharpinet et al.l (l2008 bl). It is a double optimization scheme 
aimed at isolating the model(s) that best fit the pulsation periods 
identified in this star. This is done by minimizing a X^'^yW merit 
function of general form 

^obs , pi ^ pi \ 2 

where A^obs is the number of observed periods and c, a weight 
that can be associated to each pair of observed/computed periods 
for a given model. The minimization is done both at 
the level of matching the periods and in model parameter space. 
We refer the reader to lCharpinet et al.l(i2008bl l for further de tails. 
The main difference resides in the fact that the stellar structures 
are now calculated using 3G models while in the previous study 
2G envelope structures have been used for the analysis. This has 
two consequences. First, the natural parameters needed to spec- 
ify a 3G model (and therefore the parameter space to explore) 
are different. As indicated above, these parameters are the total 



stellar mass M*, the fractional mass of the outer hydrogen-rich 
envelope log(Menv/A^*) also referred to as \ogq{H), the frac- 
tional mass of the convective core log(Mcore/A/*) or log^(core), 
and the chemical composition in the core (under the constraint 
X(He) +X{C) +X{0) = 1). Second, the external constraints 
provided by spectroscopy (Teff and \ogg) can no longer be used 
to define a priori the relevant range for the explored parame- 
ter space. This is because Tgff and \ogg are no longer natural 
parameters of the 3G stellar models, contrary to 2G envelope 
structures. Instead, these quantities now depend on the 4 param- 
eters mentioned above and their values are known only a poste- 
riori, after a model is converged. Since spectroscopic constraints 
are essential to guide the search for a meaningful asteroseismic 
solution in the vast parameter space and avoid the multiplica- 
tion of inconsistent solutions, the adopted solution is to incorpo- 
rate these constraints within the optimization procedure itself, by 
eliminating de facto during the search the models in parameter 
space that differ too much in T^^f and/or \ogg from the spectro- 
scopic values (by applying a correction factor to the 5^ value). In 
what follows we used the updated atmospheric values adopted in 
the previous section with a tolerance of 3(7 (see Table [T]i. This 
approach ensures by construction a fair consistency with spec- 
troscopy but of course there is no guarantee, a priori, that a good 
period fit exists within these constraints. 



The search for best-fit solutions was launched in the largest 
possible parameter space relevant for sdB stars: 0.30 < M^/M© 

< 0.70, -5.0 < log^(H) EE log(Me„v/M,) < -2.0, -0.40 < 
log^(core) = log(l -M^om/M^) < -0.15, and < X,ore(C+0) 

< 0.99, where Xcore(CH-O) is the fractional part of carbon and 
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0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.43 0.44 0.45 0.46 0.47 0.48 0.49 0.50 0.51 

Mass Mass 

Fig. 4. Left panel: map of the projected merit function (on a log scale) along the M*- log q{Y\) plane. Right panel: zoom of the left 
panel in the region of interest. At each M*, log ^(H) position, the value given is the projected merit function, i.e., the minimum of 
log that can be found among the values obtained for all log <7(core) and Xcore(C+0). The best-fit model is indicated by an asterisk 
on the figures, while the model corresponding to the "second" best-fit model is indicated by a cross. White contours show regions 
where the frequency fits have values within, respectively, the la, 2(7, and 3(7 confidence levels relative to the best-fit solution. 



oxygen in the convective coreQ. The constraints on log g'(H) 
and M* rely on expectations from stellar mod eling and various 
formation scenarios for hot B subdwarfs (see Han et al. "2002', 
[2003), whereas the range for the core size is loosely inspired by 
horiz ontal branch stellar evolution calculations (IDorman et al.l 
Il993h . 

We assume that the PG 1336—018 system has reached full 
spin-orbit synchronism, such that the sdB component rotates as 
a soUd body wi th P,oi = Poxh = lAlAid, h (Kilkenny et al. 2000; 
lKilkennvll201 Ih . This hypothesis relies o n very firm grounds, as 
demonstrated bv I Charpinet et al.) (l2008bl) . The rotationally split 
pulsation modes are calculated within the first-order perturbative 
approach. The importance of higher-order perturbation effects 
and tidal deform ation due to the M dwarf companion have been 
also discussed in lCharpinet et al.l (l2008bh . showing that they do 
not alter in a significant way the asteroseismic analysis, at the 
precision of the current modeling. 

As in ICharpinet et al. (2008 bi), we use 25 o ut of the 28 pe- 
riods listed in Table 4 of .Kilkennv et alj (I2003D . leaving aside 
3 periods (/lo, /n, and /g) that may be spurious. All theoretical 
modes of degree ^ = 0, 1,2 and 4 are considered in the 90—230 s 
period range. Those with ^ = 3 are explicitly excluded, having an 
extremely low visibility in the optical domain due to cancellation 
effect s because of their specific surface geometry (i Randall et alj 
I2OOI . As an additional constraint, we also limited the associ- 
ation of the dominant modes f\ to fs to I < 2, while keeping 
it open for the lower amplitude modes. This assumption sim- 
ply filters out eventual, unconvincing solutions where some ob- 



We have found that theoretical periods are not much sensitive to the 
exact core composition of C and O. Grouping (C+O) in one parameter 
facilitates and speeds up the optimization procedure. This is equivalent 
of choosing the mass fraction of helium left in the core as the free pa- 
rameter, since X^ok (He) = 1 — Xcoie (C + O). 



served periods of highest amplitude would be associated to £ = 4 
modes. This is fully justified in view of independent constraints 
on mode identification obtained from multicolor photometry for 
some bright sdB stars, from which all the dominant modes so 
far have been identified to ^ = 0, 1 and/or 2 modes (see, e.g., 
iTremblav et ani2006H\^ Grootel et al.ll2008l) . 

Within the search domain specified, the multimodal opti- 
mization code (see Charpinet et al. 2008b for details) spotted a 
region of parameter space where, strictly speaking, two close 
but distinct families of models corresponding to minima of the 
merit function exist. The absolute minimum (S^ - 4.81) has 

= 0.4699 M0, log^(H) = -3.899, log^(core) = -0.384 
and Xcore(C + 0) = 0.609. The second slightly shallower min- 
imum (S^ ^ 5.04) has quite similar stellar parameters, with 

= 0.4795 Mq, log^(H) = -3.831, log^(core) = -0.239 and 
^core(C + 0) = 0.544. With effective temperatures of 33 038 K 
(best-fit model) and 33 053 K (second-minimum model), and 
surface gravities of 5.7713 and 5.7772, respectively. Despite this 
distinction. Fig. |4] and Fig. |5] show however that the two above 
mentioned families of solutions cannot be statistically discon- 
nected. These maps are projections of the merit function (renor- 
malized to one at its minimum and shown on logarithmic scale) 
on the 2-dimensional plane, i.e., at a given point of a map, the 
value is the minimum among all the values when the two other 
parameters are varied independently. The regions of the best-fit 
model and the second-minimum model are, indeed, connected 
at the 2(7 level (white dotted contours). It is therefore appro- 
priate to treat them as a unique family, adopting the absolute 
minimum as the model of reference that best match the peri- 
ods observed in PG 1336-018. More information can be given 
on this best-fit model, including the details of the period match 
and mode identification, but since this is not essential for our 
purposes, we relegate this information to Appendix A. The only 
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-0.40 -0.35 -0.30 -0.25 -0.20 -0.15 

log q(mc) 



Fig. 5. Map of the projected merit function (in a log scale) 
along the log <7(core)- Xcoie(C+0) plane. At each log <7(core), 
^core(C+0) position, the value given is the minimum of log 
that can be found among the values obtained for all and log 
^(H). The best-fit model is indicated by an asterisk on the fig- 
ures, while the model corresponding to the second minimum of 
the merit function is indicated by a cross. White contours show 
regions where the frequency fits have values within, respec- 
tively, the la, 2a, and 3(7 confidence levels relative to the best- 
fit solution. 



relevant point for the discussion is that the relative average dis- 
persion achieved for this best fit model is AX /X — 0.18% (where 
X i s either the period P or the frequency v), corresponding to 
AP — 0.30 s and Av = 1 1 .4 /xHz. The average precision in si- 
multaneously reproducing the observed periods of PG 1336- 
018 with the 3G models is very close to that achieved with the 
2G models (AxTx = .17%, AP = 0.27 s and Av = 10.5 nHz; 
ICharpinet et al.ll2008bl) . Hence the optimal seismic 3G model is 
essentially as good as its 2G counterpart to that respect. 

In Fig. m the solid and dotted vertical lines correspond to 
the three orbital solutions and their associated I g uncertainties 
for t he mass of the sdB component derived by IVuckovic et alj 
(12007): 0.389 ± 0.005 Mq, 0.466 ± 0.006 Mq, or 0.530 ± 
0.007 Mq. The right panel of Fig. |4]is a close-up view of the 

— logq'(H) plane in the region of interest. The valley of low 
5^ -values cle arly points toward the second orbital solution of 
IVuckovic et a l. (2007), as it was alre ady found with the 2G en- 
velope models (Charpinet et al. 2008b). Fig. |5] shows an elon- 
gated valley along the log<7(core) parameter between —0.40 and 
—0.20, which reflects the fact that the p-modes observed in PG 
1336—018 are not sensitive to the stellar innermost layers, and 
therefore the size of the core cannot be constrained. The core 
composition is however indirectly constrained through its impact 
on the overall hydrostatic equilibrium of the star which affects 
the pulsation periods and the atmospheric parameters (Teff and 
log^) for a given stellar mass and envelope mass (see evolution- 
ary tracks of sdB stars from, e.g., Dorman et al. 19931 

More quantitative statements can be made by statistically 
estimating the value for each parameter of PG 1336-018 and 



its associated error from the asteroseismic fit. We adopt here a 
new procedure differing from our former methods to estimate 
erro rs based on contour m aps and significance levels (see, 
e.g. JCharpinet et aDl2008b|) . We however note that the two ap- 
proaches are essentially equivalent and provide quantitatively 
comparable error estimates (see below). We calculate the like- 
lihood function 

— 1 c2 

^{ai,a2,a^,a4) °^ e (2) 

from the J^-type merit function, S^{ai,a2,a$,a4), that has been 
sampled by the optimization code during the search for the best- 
fit models along with additional grid calculations covering the 
regions of interest in parameter space. For each parameter of in- 
terest, say fli, this function is integrated over the full parameter 
range covered by the other free parameters, thus defining a den- 
sity of probability function for the chosen parameter: 

^(ai)dfli 0= dfli y^^(fli,fl2,fl3,a4)dfl2da3dfl4 . (3) 

This density of probability function is then normalized assuming 
that the probability is equal to 1 that the value of oi is in the 
range specified for the search of a solution. In other words, the 
normalization factor is such that 

J ^(fli)dfli - 1 (4) 

over the allowed parameter range. This method permits the con- 
struction of histograms for the probability distribution of each 
primary model parameter Secondary model parameters, like 
Teff, logg, the radius R, and the luminosity L can also be eval- 
uated in a similar way. Let bi be one of these secondary param- 
eters. The density probability function is then 

[^{b\)Abi °^ Abi jjjj ^(fli,fl2,fl3,fl4;^)dflidfl2da3dfl4 (5) 

where the integration is done with the additional constraint that 
for a given set of primary parameters, the corresponding value 
b of the secondary parameter considered must be within bi and 
bi + dbi . Again, the density of probability for the secondary pa- 
rameters is normalized such that 

j .^{bi)dbi = 1 . (6) 

Figure |6] shows the histogram obtained for the probability 
distribution of the total mass of PG 1336-018 inferred from as- 
teroseismology following the procedure described above. The 
range shown in this diagram only covers the region of interest, 
where a peak in the distribution is found (there is nothing below 
0.37 Mq and above 0.56 Mo) and where the 3 orbital solutions 
for the mass of the sdB component derived by IVuckovic et al.l 
(I2007h can be represented (filled circles with a horizontal dot- 
ted line indicating the la error). Clearly, we find a very striking 
correspondence between the mass value derived from asteroseis- 
mology based on the new 3G models and the second orbital so- 
lution from Vuckovic et al. (2007). The two other orbital models 
giving respectively a much lower mass or a much higher mass 
for the sdB component are totally inconsistent with the seis- 
mic models and can be clearly excluded. From the histogram 
of the probability distribution, we can draw a statistical infer- 
ence for the mass of the sdB star. The region between the two 
vertical lines (in red) contains 68.3% of the probability distri- 
bution and therefore defines, by convention, the la range for 
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Fig. 6. Probability density function for the stellar mass from as- 
teroseismology. The filled circles with the dotted lines are the 
3 orbital solutions for the mass of the sdB component proposed 
by IVuckovic et al.l (l2007h with their la uncertainties. The red 
hatched part between the two vertical solid red lines defines the 
Icr range, containing 68.3% of the mass distribution. The blue 
vertical dashed line indicates the mass of the optimal model so- 
lution of lowest 5^-value. 



Fig. 7. Probability density function for the surface gravity log^ 
from asteroseismology. The filled circle with the solid line and 
dotted line indicate the value derived from spectroscopy with its 
1(7 and 3(7 uncertainties, respectively. The red hatched region 
between the two vertical solid red lines is the 1 a range, contain- 
ing 68.3% of the log^ distribution. The blue vertical dashed line 
corresponds to the log^ value of the optimal model solution of 
lowest 5^- value. 



the value of the represented parameter This range covers from 
0.465 to 0.477 Mq, or equivalently we infer that, statistically, 

==0.471 ±0.006 Me for the sdB star in PG 1336-018, taking 
the central value of the 1 a range as the referenc e for simplic- 
ity. We recall that the relevant orbital solution of IVuckovic et alj 
(120071) yiel ded M^, = 0.466 ±0.006 and that the seismic solution 
obtained bv lCharpinet et al.l (l2008bl) based on 2G models led to 

= 0.459 ±0.005. Interestingly, the mean and median values 
of the mass distribution are equal to 0.472 and 0.473 Mq, 
respectively. The small differences in these values indicate that 
the distribution is close to, but not exactly, a normal (gaussian) 
distribution, as can be easily seen in Fig. |6] The small visible 
"distortion" is to be related to the shape of the function illus- 
trated in the right panel of Fig. |4]for instance. Finally, we note 
that the best-fit seismic model (of minimum S'^-value) shown as 
a vertical blue dashed line lies within the 1 a range of the derived 
statistical mass distribution, meaning that this optimal model is 
not an outlier of the distribution and is representative of the star, 
as far as this parameter is concerned. 

In the spirit of confronting the results of asteroseismology 
with other techniques, two additional parameters are of specific 
interest. First, the probability distribution obtained for the sur- 
face gravity \ogg (Fig.lT) can be compared with both the val- 
ues derived from orbital modeling and spectroscopy. This dis- 
tribution shows a very narrow peak with a 1(7 interval (con- 
taining 68.3% of the distribution) between 5.768 and 5.781, or 
a statistically inferred value of \ogg = 5.775 ±0.007 (taking 
the centre of the 68.3% significance interval as the reference 
value). The comparison with other determinations shows a per- 



fect agreement, within the 1 a errors of each technique, although 
the value determined from asteroseismology is much more pre- 
cise. In particular we recall that the relevant orbital solution 
gives \ogg = 5.77 ± 0.06 (filled triangle and horizontal dashed 
line in Fig. |7]l and the spectroscopic analysis (Section 2) leads 
to \ogg = 5.771 ±0.015 (filled circle and horizontal solid line 
in Fig. IT). The former asteroseismi c solution obtained wi th 2G 
models gave log^ = 5.739 ± 0.002 dCharpinet et al.ll2008bi) . also 
consistent with the other techniques (at the achieved precision), 
but somewhat shifted compared to the new estimate based on 3G 
models. Of course, the spectroscopic value was used to constrain 
the search of an asteroseismic solution. However, Fig. |7] shows 
that the adopted 3(7 tolerance for this constraint (dotted hori- 
zontal line) is much larger than the obtained probability distri- 
bution, indicating that the search was not overly constrained by 
spectroscopy. We also checked that relaxing further the spectro- 
scopic constraints does not change at all the location and shape 
of this distribution. Finally, we note that the mean and median 
values of the distribution (5.775 and 5.776, respectively) are es- 
sentially equal and the derived asteroseismic values for logg are 
normally distributed to a very good approximation. Moreover, 
the log^ value for the optimal model (the vertical dashed line in 
blue) is representative of the distribution, as it is found within 
the 1(7 range. 

The other interesting parameter for comparison purposes is 
the radius R of the star (Fig. |8]l derived from both asteroseis- 
mology and the analysis of the orbital light curve of the system 
(in particular the eclipses). The distribution for R also shows a 
narrow peak with a 1(7 interval (containing 68.3% of the dis- 
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Fig. 8. Probability density function for the stellar radius from 
asteroseismology. The filled circle with the solid horizontal line 
indic ates the 1 ff er ror associated to the relevant orbital solution 
from lVuckovic et a l. (2007). The red hatched region between the 
two vertical solid red lines is the la range containing 68.3% of 
the radius distribution. The blue vertical dashed line corresponds 
to the radius of the optimal model solution of lowest 5^-value. 



tribution) between 0.1465 Rq and 0.1483 7?©, or a statistically 
inferred value of R = 0.1474 ±0.0009 (taking the centre of the 
68.3% significance interval as the reference value). This is a very 
precise estimate which is in total agreement with the value de- 
rived for the relevant orbital solution of Vuckovic et al. (2007). 
The latter is however an order of magnitude less precise with an 
estimated value of 7? = 0. 15 ± 0.01 Rq. The previous asteroseis- 
mic solution obtained w ith 2G models led to /? = 0.151 ±0.001 
dCharpinet et al.l l2008b). also consistent with the orbital value, 
but shifted notably (in view of to the high precision measure- 
ments) compared to the new estimate based on 3G models. 
Again, we note that the mean and median values of the distribu- 
tion (0.1474 and 0.1475, respectively) are nearly equal and the 
statistically derived asteroseismic value for R essentially follows 
a normal distribution. We also note that radius of the optimal 
model (the vertical dashed line in blue) is representative of this 
distribution (as it is within the 1(7 range). 

Asteroseismology provides much more information on the 
star than the three parameters discussed above. For complete- 
ness, we show in Fig. |9]the histograms obtained for the prob- 
ability distributions of the other derived parameters. In a nut- 
shell. Fig. |9] shows that several parameters, namely the mass of 
the H-rich envelope logq{H), the mass fraction of C+O in the 
core Xcore(CH-O) (or equivalently the mass fraction of helium left 
there), and the luminosity of the star are rather well measured by 
asteroseismology at a quite interesting precision (see Table|2]for 
the corresponding values and error estimates). It is important to 
recall that at least two of the above mentioned parameters cannot 
be measured by other means. For the three parameters discussed 
above, we also point out that the histograms shown in Fig. |9]are 



all close to normal (gaussian) distributions. In contrast, the size 
of the mixed core log^(core) remains unconstrained, as the cor- 
responding panel in Fig. |9] clearly shows. The distribution for 
this parameter is almost uniform and no useful information can 
be obtained. This is not surprising since the p-modes observed 
in PG 1336-018 are essentially envelope modes that are largely 
insensitive to the detailed structure of the deepest regions, in- 
cluding in particular the core. It is only with g-mode pulsators 
that this quantity can be effectively measured. The last param- 
eter is the effective temperature Teff which is well know n to be 
poo rly measured by the p-modes in sdB stars (Charpin et et al.l 
|2005). Indeed, in that case the corresponding panel in Fig. |9] 
shows that it is mainly the spectroscopic constraint imposed dur- 
ing the search for an asteroseismic solution that limits the range 
of the distribution. Finally, we stress that all the parameter val- 
ues characterizing the optimal model that provides the absolute 
best match to the pulsation periods of PG 1336-018 fall within 
(or very close to) the 1 a range of their corresponding statisti- 
cal distributions (Fig.|6]|7][8] and|9l blue vertical dashed line). It 
proves that this optimal model is not an outlier of the distribution 
of models that could potentially fit reasonably well the pulsation 
periods, but instead a realization of the statistics which is repre- 
sentative of the hot B subdwarf in PG 1336-018 and can safely 
be considered as the seismic reference for this star. 

3.4. Comparison between various estimations 

We summarize in Table |2] all the values inferred for the param- 
eters of PG 1336-018 from the various techniques employed. 
This sets the stage for a comparative discussion of the vari- 
ous results and an assessment of eventual systematic effects. 
As already mentioned, the best tested parameters are M, R, and 
logg because they have been measured quite precisely from both 
asteroseismology and orbital modeling, as well as from spec- 
troscopy for the surface gravity. For these 3 parameters, all mea- 
surements indicate a particularly strong (within the 1 a errors) 
consistency. Notably, the stellar mass is precisely and accurately 
determined from asteroseismology using the 3G models accord- 
ing to model II of Vuckovic et al. ( 20071). This consistency was 
already proven for the 2G models (iCharpinet et al.ll2008bl) . in- 
dicating that no important bias exists in the determination of 
this parameter from asteroseismology. This constitutes a very 
strong test demonstrating the reliability of the models (both 2G 
and 3G) for this type of seismic analyses. This is also a crucial 
point for the project of building an empirical mass distri bution 
of sdB stars from asteroseismology (iFontaine et al.ll20T2b . Such 
empirical mass distributions are essential in order to clarify the 
question of the formation of sdB stars, and unbiased measure- 
ments are of course a key factor 

A closer comparison of the two sets of seismic parameters 
derived from 2G and 3G models, respectively, shows that in 
practical terms no important difference exists between the two 
evaluations. The intrinsic precision of these measurements is 
very high by usual standards, but in general the differences be- 
tween the two analyses remain below the 3(7 errors (although 
several parameters can differ by more than 1 a). There is one ex- 
ception, log(Menv/A^*), showing a significant systematic effect, 
as this parameter is estimated to be somewhat larger with the 
3G models. This is related to the presence of a realistic stel- 
lar core in 3G models, which has a small but noticeable im- 
pact on the frequencies of p-modes, especially for degrees I — 
(Charpinet et al. 2002). These small frequency shifts slightly 
modify mode trapping, which depends on the positions of the 
eigenfunctions nodes relative to the envelope He/H transition 
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Fig. 9. Probability density functions derived from asteroseismology for the parameters loq ^(H), log ^(core), Xcore(C+0), T^ff, and 
L*/L0. The red hatched regions between the two vertical solid red lines shows the IcT range containing 68.3% of the distribution. 
The blue vertical dashed line corresponds to the value of the optimal model. In the panel showing the distribution for Tiff, the filled 
circles is the spectroscopic value with its 1 a (solid horizontal line) and 3(7 (dotted horizontal line) uncertainties. 
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Table 2. Structural parameters of PG 1336—018 derived from asteroseismology, spectroscopy, and orbital light curve analysis. 



Quantity 


Asteroseismology 


Spectroscopy 


Orbital light curve modeling 




Ihis study 


L2008 ' 


Misc. 


Vuckovic et al. (2007) 




(3G models) 


(2G models) 


(see Sect.im 


Model I 


Model 11 


Model III 


R/Rq 
logg 


0.471 ±0.006 (1.3%) 
0.1474 ±0.0009 (0.6%) 
5.775 ± 0.007 (0.1%) 


0.459 ± 0.005 
0.151 ±0.001 
5.739 ± 0.002 


5.771 ± 0.015 


0.389 ± 0.005 
0.14 ±0.01 
5.74 ± 0.05 


0.466 ± 0.006 
0.15 ±0.01 
5.77 ± 0.06 


0.530 ± 0.007 
0.15 ±0.01 
5.79 ± 0.07 


Teff (K) 

logA'(He)/A'(H) 


32850 ±175 (0.5%) 


32 740 ± 400 


32 807 ± 82 
-2.918 ±0.089 








log(Menv/M,) 

log(l-MeoreM) 

^core(C+0) 

log(L/Le) 


-3.83±0.06(1.6%) 
unconstrained 
0.58 ±0.06 (10%) 
22.9 ±0.6 (2.6%) 


-4.54 ± 0.07 
23.3±1.5 










My(g,reff,M,) 

d{V,Mv) (pc) 

/"rot (S)* 

Veq(/'iot,i?) (km/s) 

in 

V sin; (km/s) 


4.60 ± 0.04 (0.8%) 
571* ±35 (6.1%) 

8727.7823 ±0.0002 
73.9 ±0.5 (0.6%) 

72.9 ±0.5 (0.6%) 


4.49 ±0.04 
619±38 

75.9±0.6 

74.9 ±0.6 


<79* 


80.67 ±0.06 


80.67 ±0.06 


80.67 ±0.06 



^ From lCharpinet et al.lll2008bh . 

* Corrected for a reddening of E{B — V)= 0.021 ± 0.009. The 2008 distance estimate was not corrected for reddening. 

* Sync hronized at the orb ital period value tfrom lKTlkennvll^oTTh . 

* From lGeier eT^ (120101) . 



dCharpinet et al.ll2000l) . Mode trapping therefore primarily de- 
termines the value of log q(H), which thus slightly but signifi- 
cantly changes from 2G to 3G models. Overall, if the improved 
3G complete stellar models of sdB stars are now to be preferred 
for asteroseismology purposes, this comparison underlines that 
the former 2G envelope models already provided quite robust 
estimations of the structural parameters of PG1336-018. By ex- 
tension, this considerably increases our confidence in the seismic 
parameters derived for all the other p-mode sdB pulsators based 
on these 2G models. 

Finally, we point out that this reliability check for the 3G 
models using PG1336-018 as a privileged testbed remains par- 
tial in the sense that only p-modes are observed in this pulsat- 
ing sdB star. An extension of this test to a g-mode pulsator in a 
close eclipsing binary may become possible with the sdB-ndM 
system 2M 1938+4603 recently disc overed by the JCepJer space 
mission (0ste nsen et al . 2010). B arlo wet al.l (l20 1 2h proposed an 
orbital solution from Kepler eclipse timings, assuming a perfect 
circular orbit. The mass ratio, and therefore the masses of the 
two components, are however very sensitive to the eccentricity. 
The light curve modeling of 2M 1938+4603 still has to be car- 
ried out in a fully consist ent way, i.e., by using additional mul- 
ticolor photometry like in I Vuckovic et al.l (l2007l) . On the astero- 
seismology front, the sdB component of 2M 1938+4603 exhibits 
an extremely rich pulsation spectrum of both p- and g-modes. 
Moreover, the star is presumably a fairly fast rotator if synchro- 
nization is achieved (which is likely considering the very short 
orbital period of 3.024 hours). For ^-modes that have periods 
comparable to the expected rotation period of the sdB star (if 
synchronized), the simple perturbative approach to compute the 
pulsation spectrum can no longer be applied. One must then in- 
vestigate the problem with more sophisticated methods based 
either on the traditional approximation or on nonperturbative 
treatments to compute the effects of rotation on stellar pulsations 



(iReese et al.l2009HOuazzani et alJ20T2h . This star therefore cur- 
rently poses a difficult challenge for the seismic modeling, but 
could become exploitable in the future, as work being done in 
that sense progresses. 

3.5. Impact of known model uncertainties on the seismic 
inference of stellar parameters 

The static 2G and 3G models developed for detailed asteroseis- 
mic studies of sdB stars include a simplified treatment of the 
effects of diffusion. Diffusion, along with other potentially com- 
peting mixing mechanisms, is a process that can modify the 
chemical stratification inside the star over time. It is therefore 
more difficult to implement in a strategy based on static stellar 
structures and approximate ways of dealing with it must be used. 
Diffusion will most notably have two main effects. First, it will 
control the heavy metal abundance profiles in the envelope of a 
sdB star under the pressure of selective radiative forces (radiative 
levitation). Second, it will also modify (smooth), as time passes, 
the profile of the chemical transition at the interface between 
the He mantle and the H-rich envelope. We investigate below 
the impact of uncertainties associated with these two effects on 
the seismic determination of stellar parameters. We also investi- 
gate the impact of uncertainties on He-burning rates, namely, the 
triple-a and '^C(a,7)^^0 nuclear reactions. 

3.5.1. Impact of the iron composition profile 

In the simplified approach used in our models, only iron is con- 
sidered (as the main contributor to the Rosseland mean opacity) 
and the stratification adopted is assumed to be the profile com- 
puted at diffusive equilibrium, i.e., equilibrium between radia- 
tive levitation and gravitational settling. This approximation is 
both practical (in a hydrostatic modeling context) and physically 
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justified by the fact that diffusion timescales are short in the rel- 
evant region of the envelope, thus presumably leading almost 
"instantly" (relati ve to the secu lar evolut ion timescale) to a n 
equilibrium state (iFontaine et al. 2006; C harpinet et al.ll2008ah . 
Of course, in reality, diffusive equilibrium is never completely 
reached and other competing phenomena can potentially affect 
the metal distribution such as stellar winds, turbulent mixing 
(induced, e.g., by differential rotation), and/or double diffusive 
(thermohaline) convection due to /i-gradient inversions. These 
effects add to the difficulty of incorporating diffusion in static 
models and are quite uncertain in their details in evolutionary 
calculations as well. Therefore, strong uncertainties remain on 
the exact abundance distributions of chemical species inside the 
envelope of sdB stars. Because heavy elements (in particular 
iron) are the dominant source for the gas opacity, their distri- 
bution will affect the thermal structure of the star and conse- 
quently modify the periods of oscillation modes. However, this 
assessment has to be made more quantitative in order to truly 
appreciate the impact in the context of deriving structural pa- 
rameters by asteroseismology. We demonstrate below that the 
uncertainty associated with the composition profiles built up by 
diffusion and other processes in competition is not presently a 
strong issue when it comes to infer structural parameters of sdB 
stars from their pulsation spectrum. 
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Fig. 10. Upper panel: Iron abundance profiles used in four, oth- 
erwise identical, test models. The black, red, blue, and green 
curves correspond to the profile at diffusive equilibrium between 
gravitational settling and radiative levitation, the profile when 
the amount of levitating iron is decreased by a factor of two, the 
profile after a decrease of factor of four, and the profile assum- 
ing a uniform iron distribution in solar proportion. Lower panel: 
Corresponding profiles for the Brunt- Vaisala (solid curves) and 
Lamb (dotted curves for I — 1 and I — 2) frequencies. 



The first experiment is based on a set of representative 3G 
models built with the same input parameters but assuming dif- 



ferent kon abundance profiles (see Fig.fTOt. The first model rep- 
resents our standard 3G structures that use iron profiles expected 
at diffusive equilibrium. Two more models are constructed by 
artificially reducing the amount of levitating iron by a factor of 
two and four, respectively. A fourth model is built with a uniform 
iron abundance distribution in solar proportion. The purpose of 
this series of models is not to explore exhaustively all imagin- 
able abundance profiles, which would of course be impossible 
to achieve. Instead, it allows us to sample a range of situations 
between two extreme cases where, on one side, no diffusion oc- 
curs at all (the uniform solar model) and, on the other side, full 
unperturbed diffusive equilibrium is reached (the standard model 
at diffusive equilibrium). We expect that the true profile is some- 
where in between these two situations, the following results thus 
providing an accurate view of the magnitude of the induced ef- 
fects in the context of asteroseismology. 

Clearly, from the noticeable differences observed in the pro- 
files of the Brunt-Vaisala and Lamb frequencies in the corre- 
sponding models (lower panel of FigfTOli. the different iron pro- 
files should have some impact on the oscillation modes. Table|3] 
confirms this by providing the adiabatic periods of the low- 
degree (^ = — 4) p-modes computed for each of the four mod- 
els illustrated in Fig. [10] The changes in the iron profile generate 
noticeable differences in periods between the standard model at 
diffusive equilibrium (dubbed "1 • A^(Fe)/A^(H)" in TableEJ and 
the other models. These however remain relatively limited for 
models were the amount of levitating iron has been decreased by 
a factor of two and four, suggesting that the impact of the exact 
chemical composition in the sampled range should remain small. 
Not unexpectedly, the model with a uniform Fe profile in solar 
proportion shows larger variations, possibly leading to more sig- 
nificant systematic drifts in the model parameters derived from 
asteroseismology. 

In order to quantify more precisely the real impact of these 
period changes, we reanalyzed the star PG 1336-018 follow- 
ing the same procedure as before (see Sect. 13. 3K but using the 
various model assumptions about iron composition profiles dis- 
cussed above. The results obtained with the standard 3G mod- 
els (using an iron abundance profile at diffusive equilibrium), as 
discussed in Section [331 are summarized for convenience in the 
fifth column of Table H) The probability density functions ob- 
tained for the main structural parameters of the star depending 
on the abundance profiles assumed for iron in the stellar models 
are shown in Fig. IB. II IB. 21 and IB.3l of Appendix iBl The corre- 
sponding statistical inferences for these parameters are given in 
Table |4] For each parameter, we also indicate in this table the 
largest systematic drift relative to the standard model (obtained 
by comparing with the model having a uniform and solar iron 
abundance) and the drift which is, in our view, more likely rep- 
resentative of the uncertainties associated with the metal abun- 
dance profiles in the envelope of sdB stars. 

The main and obvious result of this exercise is that the over- 
all impact of changing substantially the iron abundance profile is 
very small when it comes to infer the structural parameters of the 
star from asteroseismology. Systematic effects in the derived val- 
ues, when they exist, remain quite limited in amplitude relative 
to the current precision of the measurements. The most affected 
parameters are: the radius with a drift of —0.0068 ± 0.0014 Rq 
(4.9(7) in the worst case, but more reasonably a systematic ef- 
fect of -0.0031 ±0.0015 Rq (2.1a); the surface gravity \ogg 
with a drift of -1-0.039 ± 0.010 dex (3.9ff) in the worst case, and 
more reasonably +0.020 ±0.011 dex (1.8(7); the H-rich enve- 
lope mass \ogq{H) with a drift of -0.28 ± 0.010 dex (2.8(7) in 
the worst case, but more reasonably no significant effect; and the 
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Table 3. Periods computed for a series of test models with various iron abundance or He/H transition profiles, and modified nuclear 
reaction rates. 



Fe prof." : 
He/H prof.' 



jV(Fe) t 

■IvpIT 
sharp^ 



1 N(Fe) 

sharp 



4 ■ ivtHy 

sharp 



Solar Fe 
sharp 



1 A^(Fc) 

smooth 



N{Fc) 

sharp 



Nucl. rates'" : 


CFoS ' 


CF88 


CF88 


CF88 


CF88 




i 


k 






Period (s) 









5 


95.361 


93.021 


Q1 OQI 


QS QQ7 
yo.yy 1 


95.558 


97.282 




A 


110.383 


107.854 


105.720 


114.996 


110.289 


112.633 





3 


121.878 


120.793 


119.866 


125.564 


122.073 


124.161 





2 


143.687 


139.815 


136.578 


152.200 


143.704 


146.740 





1 


167.983 


167.433 


167.018 


169.869 


167.967 


170.929 








181.622 


180.359 


179.550 


215.586 


181.500 


185.170 


} 


6 


91.986 


89.748 


o / .uu / 




91.911 


93.815 




5 


101.879 


100.347 


99.033 


104.349 


102.002 


103.844 


1 


4 


114.405 


111.923 


110.009 


120.653 


114.567 


116.759 


1 


3 


138.376 


135.678 


132.900 


142.634 


138.265 


141.102 


1 


2 


149.437 


147.657 


146.628 


155.379 


149.606 


152.381 


1 


1 


180.242 


178.677 


177.624 


214.649 


180.120 


183.837 


2 


5 


96.122 


93.970 




yy .J 


96.370 


98.043 


2 


4 


111.296 


108.611 


106.408 


116.728 


111.270 


113.597 


2 


3 


125.484 


124.479 


123.462 


128.355 


125.595 


127.811 


2 


2 


144.146 


140.426 


137.428 


152.362 


144.254 


147.202 


2 


1 


176.711 


175.225 


174.182 


183.319 


176.626 


180.125 


2 





186.021 


185.753 


185.606 


213.268 


186.069 


189.270 


3 


5 


93.563 


91.148 


89.037 


97.147 


93.663 


95.449 


3 


4 


106.995 


104.923 


103.075 


110.121 


106.922 


109.112 


3 


3 


117.429 


115.788 


114.560 


122.408 


117.701 


119.732 


3 


2 


1/11 Q/IQ 


1 JO.U /O 


134.813 


149.620 


1 A 1 QMQ 
141 .VUo 


1/1/1 Q/IM 
144.o4U 


3 


1 


163.282 


162.570 


162.053 


165.622 


163.304 


166.224 


3 





180.063 


178.710 


177.833 


210.379 


180.184 


183.593 


4 


5 


91.650 


89.378 


87.294 


94.936 


91.627 


93.479 


4 


4 


102.519 


100.850 


99.374 


105.011 


102.590 


104.508 


4 


3 


1 14.444 


112.144 


110.397 


120.172 


114.697 


116.777 


4 


2 


139.680 


136.190 


133.039 


146.380 


139.725 


142.574 


4 


1 


156.689 


155.694 


155.019 


159.771 


156.786 


159.599 


4 





178.699 


177.224 


176.246 


190.230 


178.984 


182.218 



' Standard model of reference 

" Iron abundance profile (see Section l3.5.U 
Envelope He/H transition profile (see l3.5.2"l i 
Nuclear reaction rates for the triple- a reaction (see l3.5.3l l 



luminosity L with a drift of —2.3 ± 0.9 L© (2.6a) in the worst 
case, and more reasonably —1.1 ±0.9 L© (1.2a). Of particular 
interest, we find that the mass determination does not suffer of 
any systematics at all (the value +0.003 ±0.0011 is consistent 
with no drift), thus strengthening the robustness of mass determi- 
nations of sdB stars based on asteroseismology. Consequently, 
the exact stratification of heavy metals that results from the com- 
petition of various diffusive and mixing processes occurring in 
the star and which is subject to large uncertainties is clearly not 
a strong issue at the currently achieved precision for the stellar 
parameters derived from asteroseismology. 

Nevertheless, we point out that the test case of PG 1336-018 
is further enlightening here, because it provides a comparison 



with parameters independently derived from other techniques. 
Quite interestingly. Table |4] and Fig. IB. II IB. 21 and IB. 31 show 
that our standard models assuming a nonuniform iron abundance 
profile at diffusive equilibrium lead to a seismic solution that 
is the closest to the values derived for \ogg and R from spec- 
troscopy and orbital light curve analysis. Moreover, the best-fit 
solution to the observed periods is found to be noticeably better 
with these models, as illustrated in Table|4]where the optimal 
values, mean relative period (or frequency) differences AX /X, 
and mean frequency differences Av achieved for the optimal fits 
are compared. These facts suggest that the standard models are 
the most realistic structures among this batch of test models to 
represent the oscillation properties of PG 1336-018 and that the 
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Table 4. Parameters derived for PG 1336—018 using various iron abundance profiles. 



Parameter 


Uniform/Solar 


i . ^(Fe) 


1 N(Fe) 

' Ml'U\ 

^ J V ( rl j 




largest drift 


likely drift 


M/Mq 


0.471 ±0.009 


0.474 ±0.009 


0.474 ±0.009 


0.471 ±0.006 


+0.003 ±0.011 


+0.003 ±0.011 






C TOC 1 A AAO 

5.795 ±0.00o 


C TOjC 1 A AAA 


C TTC 1 A AAT 

5.775 ±0.007 


1 A AQA 1 A A1A 

+0.039±0.010 


1 A AO A 1 A Al 1 


"/«© 


n 1 ACif, -1- 00 1 1 


n 1 ziA^ + n nn 1 9 

U . i ft .3 U . UU i ^ 


1 ASQ -1-0 001 9 


fl 1 474 + f) flflflQ 


fl nn(S8 + n nnM 

U.UUUO 31 U.UUiH 


A 00^ 1 001 s 


logq{H) 


-4.11 ±0.08 


-3.70 ±0.09 


-3.74 ±0.09 


-3.83 ±0.06 


-0.28 ±0.10 


+0.07 + 0.11 


X{C + 0) 


0.46 ±0.10 


0.48 ±0.10 


0.51 ±0.09 


0.58 ±0.06 


-0.12±0.12 


-0.10±0.12 


L/Lq 


20.6 ±0.7 


21.8±0.7 


22.3 ±0.8 


22.9 ±0.6 


-2.3 ±0.9 


-1.1±0.9 


52 (opt.) 


5.95 


5.93 


5.46 


4.81 






AX/X(%) 


0.24 


0.21 


0.19 


0.18 






Av (MHz) 


14.07 


12.75 


11.43 


11.37 







iron profile at diffusive equilibrium is currently the approxima- 
tion that is the closest to the real iron distribution in the star We 
however stress that improving the treatment of diffusion in the 
next generations of sdB models in order to obtain closer fits to 
the observed periods and increase the precision of the measured 
quantities is among our long term objectives. 



3.5.2. Impact of the envelope He/H transition profile 

As a sdB star evolves, diffusion affects the profile of the chem- 
ical transition between the H-rich envelope and the He mantle. 
From a relatively sharp interface at the beginning of core He- 
burning, the transition will slowly spread over time leading to 
smoother profiles (see, e.g. jHu et al.l l201 iT). 

In our static 2G and 3G models, the profile of the He/H 
transition zone is controlled through a single parameter that 
sets the maximum slope of the composition gradient. This pa- 
rameter has been c alibrated based on evolutionary models from 
iDorman et al. I (11993.) that do not incorporate diffusion of helium. 
We therefore traditionally assume a sharp He/H transition in our 
standard stellar models used for asteroseismology and the im- 
pact of this assumption on the seismic solutions is a legitimate 
question. We demonstrate in the experiment below that the shape 
of this chemical transition indeed slightly affects the periods of 
the relevant oscillation modes, but at a level that has no practi- 
cal consequence for deriving structural parameters by asteroseis- 
mology. 

We first considered two models with the same structural pa- 
rameters (and our usual nonuniform iron abundance distribution 
from diffusive equilibrium) but with different He/H transition 
profiles (see the upper panel of Fig. fTTT i. The first model, which 
is the standard model of reference (see Section l37Jt . has a sharp 
transition. The second model is built assuming a smooth transi- 
tion. The Brunt-Vaisala and Lamb frequency profiles shown in 
the lower panel of Fig.[TT]clearly illustrate that the differences in 
model structure exclusively occur in the transition region, with 
a wider Ledoux spike for the model having a smooth transition. 
The profiles remain nearly unchanged in the other regions of the 
star It suggests that a rather limited impact on the pulsation peri- 
ods occurs, which is confirmed in Table|3]by comparing the third 
and seventh columns. Between the two models, the period differ- 
ences for the low-order p-modes of interest remain very small. 
It is therefore anticipated that the shape of He/H transition zone 
will have no significant impact on the structural parameters de- 
rived from asteroseismology. 



1 - 
0.8 
0.6 
0.4 
0.2 

- 



og q - log(l -m{r)/M) 




-5 -10 
log q log{ 1 -m(r)/M) 



Fig. 11. Upper panel: Helium abundance profiles used in two, 
otherwise identical, test models. The curves in black and red re- 
spectively show a sharp and a smooth He/H transition profile 
at the envelope lower boundary. Lower panel: CoiTesponding 
profiles for the Brunt-Vaisala (solid curves) and Lamb (dotted 
curves for / = 1 and / ~ 2) frequencies. 



We explicitly demonstrate the above assertion by reana- 
lyzing PG 1336-018 with models featuring a smooth He/H 
transition. The probability density functions obtained for the 
main structural parameters of the star are provided in Fig. B.4I 
(AppendixIBI). The corresponding statistical inferences for these 
parameters are given in Table |5]for comparisons with the values 
derived with our standard models assuming a sharp transition. 
No significant drift in the derived stellar parameters is observed 
and the achieved quality of the optimal period fit is compara- 
ble. Consequently, the effects induced by diffusion on the man- 
tle/envelope transition profile presently have no noticeable im- 
pact on deriving the structural parameters of sdB stars by astero- 
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Table 5. Parameters derived forPG 1336—018 using a sharp or 
a smooth He/H transition profile. 



Parameter 


shaip transition 


smooth transition 


drift 


M/Mo 


0.471 ±0.006 


0.476 ±0.009 


±0.005 ±0.011 


log« 


5.775 ±0.007 


5.777 ±0.007 


+0.002±0.010 


R/Rf 


0.1474 ±0.0009 


0.1476 ±0.0011 


±0.0002±0.0014 


\ogq(H) 


-3.83 ±0.06 


-3.79 ±0.08 


±0.04 ±0.10 


X{C + 0) 


0.58 ±0.06 


0.53 ±0.09 


-0.05 ±0.11 


L/Lr. 


22.9 ±0.6 


22.8±0.7 


-0.1 ±0.9 


(opt.) 


4.81 


4.91 




AX/X(%) 


0.18 


0.18 




Av (/jHz) 


11.37 


11.48 





seismology, considering the actual precision of the fits. This can 
safely be regarded as a high order effect. 

3.5.3. Impact of nuclear reaction rates 

In a last experiment, we investigated the effects of nuclear reac- 
tion rates on the inferred parameters of the star Our 3G stellar 
models incorporate the network of nuclear reactions (in partic- 
ular the triple- a and related reactions for He-buming) to com- 
pute the rate of energy production that will influence the overall 
thermal equi librium of the star. For our purposes, the code still 



relies on the ICau^hlan & Fowleil (Il988l) nuclear reaction rates. 



However, newer rates from lAngulo et al ] (11999) have since been 
published which, at the temperatures relevant to He-buming in 
sdB stars, increase by nearly 10% the triple- a rate and by a fac- 
tor of 2 the rather uncertain '^C(a,7)'^0 rate. 

In order to estimate the real impact of these changes on 
the seismic inferences of the stellar parameters, we constructed 
test models where the nuclear energy production rates are in- 
creased by 10% for the triple- a reaction and doubled for the 
^2C(a,7)'^0 reaction. These are referred to as the "2xCF88" 
models. The effect of this change on the pulsation periods is 
given in last column of Table [3] Compared to the model of refer- 
ence, a slight nearly uniform increase of the periods is observed. 
This reflects the fact that for the same input parameters, the 3G 
model with the modified nuclear reaction rates leads to a star 
with a slightly larger radius, and therefore slightly less compact. 



Table 6. Parameters derived for PG 1336—018 using models 
with modified reaction rates for the triple-a and '^C(a,7)'^0 
reaction rates. 



Parameter 


CF88 


2xCF88 


drift 


M/Mg 


0.471 ±0.006 


0.474 ±0.010 


±0.003 ±0.012 


log« 


5.775 ±0.007 


5.775 ±0.008 


0.00 ±0.0 11 


R/Re 


0.1474 ±0.0009 


0.1479±0.0010 


±0.0005 ±0.0013 


iogq(H) 


-3.83±0.06 


-3.84 ±0.09 


-0.01 ±0.11 


X(C + 0] 


0.58 ±0.06 


0.54 ±0.10 


-0.04 ±0.12 


L/Le 


22.9±0.6 


23.0±0.7 


±0.1 ±0.9 


(opt.) 


4.81 


5.12 




AX/X (%) 


0.18 


0.19 




Av(mHz) 


11.37 


11.52 





We reanalyzed PG 1336-018 with the "2xCF88" models, 
and the results are given in Fig. IB . 5 1 (Appendix IBTi showing the 
probability density functions obtained for the main structural 



parameters of the star The corresponding statistical inferences 
for these parameters are given in Table |6]for comparisons with 
the values derived with our standard model. No significant drift 
in the derived stellar parameters is observed and the quality of 
the optimal period fit is comparable. Cons equently, the impact 
of upd ati ng the nuclear reac tion rates from ICaughlan & Fowlen 
(Il988h to lAngulo et aP (H^Q) on the values derived from aster- 
oseismology is not expected to be important. This will be done 
in future versions of our static models aimed at detailed seismic 
analyses of sdB stars, but clearly this is currently not a major 
source of uncertainty for our current purposes. 

4. Summary and conclusion 

In this paper we presented new stellar models of so-called third 
generation (3G), aimed at performing detailed asteroseismic 
studies of hot B subdwarf pulsators. These models are static 
structures realistically representing the star from the surface to 
its center, contrary to the former generation of envelope mod- 
els of second generation (2G) that were missing the innermost 
regions. This improvement was a strong necessity in order to 
exploit with asteroseismology the ^-mode sdB pulsators that 
have been and are currently observed from space with CoRoT 
and Kepler. The new models are constructed assuming hydro- 
static and thermal equilibrium for a given set of model parame- 
ters, including a chemical composition of the core. They allow 
for the computation of accurate periods/frequencies for the p- 
modes and the g-modes (including possible mixed modes), while 
the former generation of envelope models was limited to the p- 
modes in terms of accuracy. 

In order to test the reliability of these new 3G models for as- 
teroseismology, we reanalyzed the pulsating sdB star PG1336- 
018 with them. This object is a close eclipsing sdB-ndM system 
where the sdB component is a p-mode pulsator This very rare 
configuration allows us to test some of the stellar parameters 
inferred from asteroseismology with the values obtained inde- 
pende ntly from the orbital light curve analysis of IVuckovic et al.l 
(I2007h . Our new seismic results show that a very strong con- 
sistency persists at that level, as it was shown to be the case 
with our former 2G mod els used in the seismic study of this star 
dCharpinet et al.ll2b08bl) . Our new analysis provides updated es- 
timates of the structural parameters of the sdB star in PG1336- 
018 based on more realistic models (see Table |2]l. The derived 
parameters for the star do not diff er much from those already 
inferred in ICharpinet et al.l (l2008bb . but more importantly they 
demonstrate that the asteroseismic inferences based on both 2G 
and 3G models are both precise and accurate, at least for the 
stellar parameters that could be tested, which includes the mass 
and radius of the star. This is an important result in view, for 
instance, of our current attempts to construct an empirical mass 
distribution for comparis on with predictions of various forma- 
tion channels of sdB stars (iFontai ne et al. 2012). We also demon- 
strated that known model uncertainties do not affect, at the cur- 
rent level of accuracy, the structural parameters derived from as- 
teroseismology. However, this is one of our long term goal to in- 
clude improved physics in stellar models, in particular the treat- 
ment of diffusion, in order to obtain closer fits to the observed 
pulsations and increase the precision of the derived quantities by 
asteroseismology. 

Our approach based on the minimization of a merit function 
has sometimes been criticized in the past, arguing that the best- 
fit seismic model is not necessary the most representative model 
of the star We show in this analysis, through an improved treat- 
ment and estimation of the statistical significance of the inferred 
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parameters, that the optimal model that best match the observed 
periods of PG1336-018 is not an outlier of the statistical distri- 
bution of potential solutions. It can therefore safely be consid- 
ered as the model that best represents the star in consideration, 
as we used to consider it in our previous analyses. 

Overall, we can conclude that the main goal of this paper, 
which was to present and demonstrate the basic validity of the 
new thkd generation models used for asteroseismology of sdB 
stars, is fully achieved. Such models have passed an important 
test showing that they are up to this task. 
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Appendix A: Period fit and mode identification 
from the optimal model of PG 1336-018 

The period match and mode identification obtained for the 
best-fit model of PG 1336-018 are given in Table \AJ\ The 
relative and absolute differences for each pair of the 25 ob- 
served/theoretical modes are provided, in period and frequency, 
AX/X (in %), AP (in sec), and Av (in ^Hz). On average, the rel- 
ative dispersion between the observed and theoretical periods is 
0.18%, which is similar to the best-fit model of second genera- 
tion uncovered by Charpinet et al. (2008b) that led to 0. 17%. 

Independent constraints on the mode identification in 
PG 1336—018 exist from time-resolved spectroscopy 
dVuckovic et al.l |2009). Their detailed line-profile analysis 
excluded the {£,\m\) = (3,3), (1,0), (2,1) or (2,0) modes for 
the 5435 jUHz pulsation (/g in Table lA.ll ). although data were 
too noisy to uniquely identify this pulsation with the radial or 
sectoral dipole/quadmpole modes. The 5435 /iHz /s pulsation is 
identified as the fundamental radial mode {k, £, m) = (0,0,0) in 
Table lA.ll which is consistent with the results of lVuckovic et alj 
2009, (see in particular their Fig. 13). 

It is also interesting to lo ok back at the three uncertain peri- 
ods of iKilkennv etaP ( 120031) that were not included in the seis- 
mic search: 184.04 s (/lo), 178.96 s (/n), and 173.59 s (/g). /lo 
and /s all have an acceptable counterpart in the best-fit model 
theoretical spectrum (see Table lA.ll . f\Q can be associated with 
the (l,k,m)— (2,0,— I) mode for a relative dispersion AX/X = 
-0.33%, and /g finds a counterpart at AX/X = -0.26% with 
the (l,k,m)~(4,Q,—2) mode. These matches do not significantly 
affect the overall quality of the fit. A Z = 3 mode (not shown 
here) is needed to account for the last uncertain period, 178.96 s 
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(/is)- We recall however that the existence of these periods need 
to be confirmed by additional observations. 

The mode identification of the best-fit model respects quite 
well, qualitatively, the expected amplitude hierarchy. On aver- 
age, one may expect that modes of higher I would have lower 
apparent amplitudes even though this argument can hardly be ap- 
plied to individual identifications, since intrinsic amplitudes are 
not known and may vary significantly from one mode to another 
According to the mode identification in Table lA.l"! the observed 
average amplitudes A/ are: Aq ~ 0. 145% (for / — modes), Ai = 
0.230%, A^= 0.183%, and A^= 0.081%. This significant drop 
for / = 4 modes corresponds well to the theoretical expectation 
that mo des with / — 4 have a much lower visibility than the I < 2 
modes (iRandall et alJl2005ll2007h . 
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Table A.l. Mode identification and details of the frequency fit obtained for the 
optimal solution. The mean relative dispersion of the fit is AX/X = 0.18% (or 
AP = 0.30 s and Av = 11.366 ^iHz; 52 = 4.812). 
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4 




-1 




6427.225 




155.59 












4 









6327.371 




158.04 












4 




+ 1 




6227.516 




160.58 












4 




+2 


6163.328 


6127.661 


162.25 


163.19 


-0.58 


-0.94 


+35.667 


0.0500 


ht 


4 




+3 




6027.806 




165.90 












4 




+4 


5916.110 


5927.951 


169.03 


168.69 


+0.20 


+0.34 


-11.842 


0.1000 


/ii 


4 





-4 




5969.956 




167.51 












4 





-3 




5857.737 




170.71 












4 





-2 




5745.518 


[173.59] 


174.05 










[/si 
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Table A.l. continued. 









Vobs 


Vth 


^obs 




AV IV 
ZXA 1 A 




Air 


Amplitude 


lu. 


/ 


k 


m 


(MHz) 


(juHz) 


(s) 


(s) 


(%) 


(s) 


(jUHz) 


(%) 




4 





-1 


5621.135 


5633.299 


177.90 


177.52 


+0.22 


+0.38 


-12.164 


0.0700 


/zo 


4 








5516.633 


5521.080 


181.27 


181.12 


+0.08 


+0.15 


-4.447 


0.1000 


/l2 


4 





+ 1 


5401.026 


5408.861 


185.15 


184.88 


+0.14 


+0.27 


-7.835 


0.0600 


/25 


4 





+2 




5296.642 




188.80 












4 





+3 




5184.423 




192.89 












4 





+4 




5072.205 




197.15 












4 




-4 


5444.251 


5450.925 


183.68 


183.46 


+0.12 


+0.22 


-6.674 


0.1700 


fl 


4 




-3 


5356.473 


5338.149 


186.69 


187.33 


-0.34 


-0.64 


+ 18.324 


0.0800 


/16 


4 




-4 


5218.934 


5225.373 


191.61 


191.37 


+0.12 


+0.24 


-6.439 


0.0600 


hi 


4 




-1 


5111.168 


5112.598 


195.65 


195.60 


+0.03 


+0.05 


-1.430 


0.0500 


hi 


4 









4999.822 




200.01 












4 




+ 1 


4885.198 


4887.046 


204.70 


204.62 


+0.04 


+0.08 


-1.848 


0.0700 


/18 


4 




+2 




4774.270 




209.46 












4 




+3 




4661.494 




214.52 












4 




+4 




4548.718 




219.84 
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Appendix B: Probability density functions of stellar 
parameters with known model uncertainties 
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Fig. B.l. Probability density functions derived from asteroseismology using models with a uniform distribution if iron in solar 
proportions. The parameters inferred for the sdB component of PG 1336-018 are, from the upper-left panel to the lower-right panel, 
the massM, logg, the radius R, loq ^(H), log ^(core), Xcore(CH-O), and L^/Lq. The red hatched regions between the two vertical solid 
red lines shows the la range containing 68.3% of the distribution. In the panel showing the mass, the filled circles with the dotted 
lines are the 3 orbital solutions for the mass of the sdB component proposed by Vuckovic et al. (2007) with their 1 a uncertainties. 
In the panel showing logg. The filled circle with the solid line and dotted line indicate the value derived from spectroscopy with 
its Id and 3(7 uncertainties, respectively. Finally, in the panel showing the radius. The filled circle with the solid horizontal ligg 
indicates the 1(7 error associated to the relevant orbital solution from Vuckovic et al. (2007). In each panel, the blue vertical dashed 
and dotted Unes show the seismic solution and its 1 (7 uncertainty obtained with our standard 3G models. 
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Fig.B.3. Same as Fig. IB. II but using models with a nonuniform iron abundance profile reduced by a factor of two relative to the 
value expected at diffusive equilibrium. 
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Mean : 0,4749 
Median: 0.4760 




Mean : 5,775 
Median; 5,777 




Mean : 0,1479 
Median: 0,1481 




Mean : -3,835 
Median: -3,830 




-4.4 -4.2 



-3,6 -3.6 -3,4 -3,2 

log q(H) 



Mean : 0,5396 
Median: 0,5600 




Mean : 22.96 
_ Median: 23,10 
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X(C+0) 



Fig.B.5. Same as Fig. lB.li but using models with modified rates for the triple-a and ^^C(a,7)^^0 nuclear reactions. 
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